Comment on "Four-component relativistic theory for NMR 
parameters: Unified formulation and numerical assessment of 
different approaches" [J. Chem. Phys. 130, 144102 (2009)] 

Radoslaw Szmytkowski* and Patrycja Stefariska 

Atomic Physics Division, Department of Atomic Physics and Luminescence, 
Faculty of Applied Physics and Mathematics, Gdansk University of Technology, 
Narutowicza 11/12, 80-233 Gdansk, Poland 

January 12, 2013 

Abstract 

In the paper commented on [J. Chem. Phys. 130 (2009) 144102], Cheng et al. derived a 
formula for the magnetic dipole shielding constant a for the Dirac one-electron atom in its 
ground state. That formula involves an infinite series of ratios of the Euler's gamma functions. 
We show that with some algebra the series may be expressed in terms of elementary functions. 
This leads to a simple closed- form expression for the shielding constant. 



In a recent paper [1], Cheng et al. have discussed a four-component relativistic theory of NMR 
parameters. In particular, they have applied an analytical calculation technique, based on the 
Sturmian expansion of the first-order Dirac-Coulomb Green function, found some years ago by 
one of us in Ref . [2] , to derive a formula for the magnetic dipole shielding constant a for the Dirac 
one-electron atom in its ground state. They have shown that a may be written in the form 

a = cr_i + cr+2, (1) 

where 
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with 

7 k = - (aZf. (4) 

Here Z is the nuclear charge, a is the Sommerfeld's fine-structure constant, while T(£) is the Euler's 
gamma function. 

The representation of (7+2 displayed in Eq. (3) looks formidable. It is the purpose of this 
comment to prove that the series in Eq. (3) may be summed to a closed form, leading to an 
extremely simple expression for er+ 2 . 
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To begin, we observe that with the aid of the well-known property Cr(C) = T(( + 1), we may 
write 

r(n + 72 - 7i + 2) = (n + 72 - 7i)[I> + 72 - 7i) + T{n + 72 - 71 + !)]■ (5) 
Plugging Eq. (5) into Eq. (3) yields 
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Since it is known [3, p. 36] that 
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where 2 i-\ is the hypergeometric function, Eq. (6) may be cast into the form 
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The following identity [3, p. 40]: 
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is due to Gauss. Applying it to the two 2 F\ functions appearing in Eq. (8), we find 
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Inserting Eqs. (10) and (11) into Eq. (8), after some rearrangements involving, among others, the 
identity 

7 2 2 = 7i 2 + 3, (12) 
we eventually arrive at the following expression for a +2 : 
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With no doubts, it looks much neater than the one in Eq. (3)! 

Insertion of Eqs. (2) and (13) into Eq. (1) leads to the following closed-form representation of 
the magnetic dipole shielding constant for the relativistic hydrogen-like atom in its ground state: 
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which is identical with the expression found earlier by Moore [4] , Pyper and Zhang [5] and Ivanov 
et al. [6] (after it is taken into account that the latter authors define a with the opposite sign). 
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